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Abstract 

We consider very constrained versions of the minimal supersymmetric extension of the 
Standard Model (VCMSSMs) which, in addition to constraining the scalar masses mo and 
gaugino masses mi/2 to be universal at some input scale, impose relations between the 
trilinear and bilinear soft supersymmetry breaking parameters Aq and Bq. These relations 
may be linear, as in simple minimal supergravity models, or nonlinear, as in the Giudice- 
Masiero mechanism for generating the Higgs-mixing /i term. We discuss the application of 
the electroweak vacuum conditions in VCMSSMs, which may be used to make a prediction 
for tan/? as a function of mo and mi/2 that is usually unique. We baseline the discussion 
of the parameter spaces allowed in VCMSSMs by updating the parameter space allowed 
in the CMSSM for fixed values of tan/3 with no relation between Aq and Bq assumed a 
priori, displaying contours of -Bo for a fixed input value of Aq, incorporating the latest 
CDF/DO measurement of rrit and the latest BNL measurement of — 2. We emphasize that 
phenomenological studies of the CMSSM are frequently not applicable to specific VCMSSMs, 
notably those based on minimal supergravity, which require ttlq = m^j2 as well as Aq = 
Bq + mo- We then display (mi/2,'"^o) planes for selected VCMSSMs, treating in a unified 
way the parameter regions where either a neutralino or the gravitino is the LSP. In particular, 
we examine in detail the allowed parameter space for the Giudice-Masiero model. 
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1 Introduction 

Supersymmetry is one of the most appealing extensions of the Standard Model (SM), for 
many reasons including the hierarchy problem, its necessity in string theory, unification 
of the SM gauge couplings, the suggestion of a light Higgs boson, the possibility that the 
astrophysical cold dark matter might be provided by the lightest supersymmetric particle 
(LSP) and (just possibly) the anomalous magnetic moment of the muon, — 2. However, 
supersymmetry is a general framework that accommodates many new degrees of freedom. 
The simplest possible realization of supersymmetry is the minimal supersymmetric extension 
of the SM (MSSM). Four types of supersymmetry-breaking parameters appear in the MSSM: 
scalar masses mo, gaugino masses mi/2, trilinear couplings A and a bilinear coupling B in 
the Higgs sector. In the MSSM alone, the number of free parameters associated with soft 
supersymmetry breaking exceeds 100, unless one assumes some degree of universality for 
the sparticles with different quantum numbers and flavours. In phenomenological studies 
of supersymmetry, the values of mo for the different sflavours are often constrained to be 
universal at some input GUT scale, as are the values of mi/2 for the different SM gauge group 
factors, and the A parameters corresponding to different SM Yukawa couplings, a framework 
often called the CMSSM. 

One may go even further, and assume some relation(s) between the parameters mo, mi/2, 
A, B and the gravitino mass 777.3/2- In particular, many very constrained versions of the 
MSSM (VCMSSMs) derive or postulate relations between the A and B parameters, which 
we parametrize as ^4 = Amo, B = Bm^. These relations may be linear: for example, generic 
minimal supergravity models predict that B = A — 1 [1, 2] as well as mo = m3/2, and the 
simplest Polonyi model [3] of supersymmetry breaking further predicts that \A\ = 3 — \/3 [4]. 
On the other hand, a prominent example of a nonlinear relation is 

« = ^. (1) 

A-3 ^ ^ 

which appears in the Giudice-Masiero mechanism [5] for generating the /i term. 

In the CMSSM, one may regard mo,mi/2 and A as independent parameters, and use 
the two electroweak vacuum conditions resulting from the specification of mz and the ratio 
of Higgs vacuum expectation values, tan/3, to fix and the pseudoscalar Higgs mass m^, 
which is equivalent to fixing B. As we show in this paper with some explicit examples, the 
value of B that results for any given choice of mo,mi/2 and A may not correspond to any 
plausible theoretical model. Conversely, in a VCMSSM where B is fixed in terms of A, one 
can use the electroweak vacuum conditions to predict tan /3 as a function of mo, mi/2 and A. 



1 



In a previous paper [6] , we demonstrated this type of prediction for a few specific VCMSSMs 
with hnear relations between A and B, including minimal supergravity, with the simplest 
Polonyi model as a special case. 

In this paper, we extend the previous discussion to include the Giudice-Masiero model. 
In this case, in addition to the relation (1) between B and A, the value of /i is in principle 
also predicted as 



mo 



A-3 

A 



V3 



(2) 



where X/Mpi is the coupling between a hidden sector superfield and the two Higgs superfields. 
The value of A is presumably not completely arbitrary: for example, one should probably 
require |A| ^ 0{1). This bound on |A| in turn imposes a range on the ratio ///mo for a given 
A. Since the value of |//| is an output quantity in our approach to VCMSSMs, one must 
check that |A| is not very large, which could in principle restrict the ranges of the input 
parameters. 

The first step in this paper is to discuss the application of the electroweak vacuum 
conditions. In principle, more than one value of tan/3 might be consistent with a given 
VCMSSM for some specific choice of (mo,mi/2) [7]. In practice, over large regions of the 
(mo, mi/2) we find only one solution for tan (3, as we explain in some detail. We also discuss 
the renormahzation of the input relation between A and S in a generic VCMSSM, including 
the relation between the input and electroweak-scale values of B and the one-loop threshold 
corrections at the electroweak scale. 

The second step is to update previous analyses of the CMSSM, including some updates 
in calculations of the supersymmetric particle spectrum, as well as the latest information on 
rrit, — 2, — > 57 and Bg — > /i^/i'. We demonstrate that the values of B required in the 
CMSSM for generic values of tan /3 and A do not fit within favoured VCMSSM frameworks, 
such as those based on minimal supergravity or the Giudice-Masiero model. 

We then discuss the (mo,mi/2) planes for some specific VCMSSMs, taking into account 
the fact that minimal supergravity models predict that mo = m3/2 before renormahzation, 
which is not necessarily the case in a generic CMSSM. This relation enables one to delineate 
the regions where the LSP is the lightest neutralino the lighter f or the gravitino G. 
We present unified descriptions of the x G LSP regions for some specific VCMSSMs, 
incorporating the constraints on decays of the next-to-lightest supersymmetric particle (NSP) 
into a gravitino LSP that are imposed by concordance between the Big-Bang nucleosynthesis 
(BBN) and cosmological microwave background (CMB) determinations of the baryon-to- 
entropy ratio [8, 9, 10, 11]. Finally, we discuss the Giudice-Masiero model in more detail. 
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finding that the imphed values of |A| in allowed regions of parameter space are generally 
^ particularly in the gravitino LSP region. 



2 Summary of Models of Supersymmetry Breaking 

As discussed in [6], we assume an N — 1 supergravity framework, interpreted as a low- 
energy effective field theory. In minimal supergravity models, the Kahler function K that 
describes the kinetic terms for the chiral supermultiplets $ = (C, 0), where the ( represent 
hidden-sector fields and the 0* observable-sector fields, has the form K — Ej|$'p. We denote 
derivatives of K with respect to the chiral superfields by = dK/dcf)^, etc. In the minimal 
supergravity case, we have K'^ — 0'* -\- W^W, Ki — (l)i + W*/W*, and {K'^Y^ — 5f, and the 
resulting scalar potential is (in units where the Planck mass is unity) 

V{(j), (j)*) = e-^^^'* [\W' + (j)'*W\^ - 3\W\^] . (3) 

It is then apparent that the soft supersymmetry-breaking scalar masses mo are universal at 
the input GUT scale, with [1] 

^0 = ^3/2, (4) 

where 7/13/2 is the gravitino mass and we assume that the trcc-lcvcl cosmological constant 
vanishes. If we further assume that the supcrpotential may be separated into pieces 

/ and g that are functions only of observable-sector fields 0* and hidden-sector fields (, 
respectively, then the soft supersymmetry-breaking trilinear terms Aq and bilinear terms Bq 
are also universal, and are related by [1] 

Bo ^ Ao- m3/2, (5) 

so that 

B ^ A-1, (6) 

which is one of the principal options we studied in [6] and discuss further below. 

The simplest model for local supersymmetry breaking in minimal supergravity [3] has 
just one additional chiral multiplet ( in addition to the observable matter fields (pi, with a 
supcrpotential that is separable in this so-called Polonyi field and the observable fields (pf 
W — fi'P) + 9{C)- It takes the simple form 

^(C) = ^(C + /3), (7) 
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where we impose — 2 — -\/3 to ensure that the cosmological constant vanishes. Assuming (5 
to be positive, and using (C) = -\/3 — 1, we have [4] the universal soft trihnear supersymmetry- 
breaking terms 

A = (3- V3)m3/2, (8) 
and universal bilinear soft supersymmetry-breaking terms 

B = (2-V^)m3/2, (9) 

whose consequences we explored in [6] and discuss further below. 

In the simplest version of the Giudice-Masiero (GM) mechanism [5], in addition to mini- 
mal supergravity kinetic terms in the observable and hidden sectors, and a separable super- 
potential W = f + one postulates a coupling 

K{ct>,C) 9 \C}H^H2, (10) 

where H12 are the two Higgs supermultiplets in the MSSM. Assuming that the cosmological 
constant vanishes, the term (10) generates a Higgs mixing term (2). This mechanism also 
yields the nonlinear relation between B and A given in (1), whose consequences we explore 
below. 

As already remarked, minimal supergravity models predict a relation (4) between mo and 
the gravitino mass, which is not necessarily true in the generic CMSSM. This relation enables 
us to dehneate the regions of VCMSSM parameter space where the LSP is a neutralino, the 
lighter f or the gravitino. The astrophysical and cosmological constraints on gravitino dark 
matter have been recently re-examined [8, 9], taking also into account the constraints on 
decays of the next-to-lightest supersymmetric particle (NSP) arising from comparing the 
BBN and CMB constraints on the baryon-to-entropy ratio [10, 11]. In our later discussions 
of VCMSSMs, we give unified treatments of the parts of (mi/2, mo) planes where the LSP is 
a neutralino, the lighter f and the gravitino. 

3 The Electroweak Vacuum in VCMSSMs 

In the general CMSSM, we start with the following set of input parameters defined at 
the GUT scale: mi/2, mo, Ao, Bq and the Higgs mixing parameter ^q. At tree level, the 
electroweak vacuum is specified by the following two conditions: 

2 _ 2(m^ + /x^-(mi + /x^)tan^/?) 

" (tan2/3-l) ' ^^^^ 

sin 2/3 = -2Bn/{ml + ml + 2n^), (12) 
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and the pseudoscalar neutral Higgs mass ttia is determined by 

m\ = mj + ml + 2i^'^, (13) 

where mi and 1712 are the soft supersymmetry-breaking masses for the two Higgs dou- 
blets at the electroweak scale. These as well as fi and B are assumed to be evaluated 
by renormalization-group equation (RGE) running from the input values. One may, alter- 
natively, solve for /x and B in terms of niz and tan/?: 

2 m?-mitan2/?-F|m|(l-tan2/3)-KAW 
^ ~ tan^ /3 - 1 + Ai?^ 

B/i = --(mi + m2 + 2/x^)sin2/3 + Ab (14) 

where we have now included the loop corrections Ab and A^-^'^^ required to relate the 
RGE values to the corresponding quantities evaluated at mz [12, 13, 14], and here mi^2 = 
'>Tii,2{'>Tiz) ^- In most treatments of the CMSSM, m.1/2, mo, Aq and tan/3 are taken as inputs, 
and the conditions (14) are used to determine /i, B and the CP-odd Higgs mass m^- 

As discussed in [6] , in a VCMSSM where B is determined in advance in terms of A, it is 
convenient to use the electroweak vacuum conditions (14) to determine tan/3 as a function of 
mo and mi/2 for some input value of A. However, since A^ depends on tan (3, and Ab depends 
on both fj, and tan /3 in a nonlinear way, it is not possible to write down an analytical solution 
for tan/3. Moreover, it was shown in [7] using an RGE-improved tree-level calculation for 
tan (3 in the minimal supergravity model that there may be up to three possible solutions 
for tan/3 for any given choices of mi/2,mo, and ^40. We remarked previously [6] that we 
typically find just one solution with a moderately low value of tan /3, that multiple solutions 
exist only for < 153 GeV, and that B always increases monotonically with tan /3 over the 
range 5 < tan/3 < 55 in our calculations. Thus, a given value of mi/2, mo, Ao and sgn{ii) 
always corresponds, in our analysis, to a definite value for tan /3. Since this is important for 
our treatment of VCMSSMs, we now illustrate this point in more detail. 

We show in Fig. 1 some examples of the necessary input values of Bo as functions of tan /3, 
for four representative choices of (mo,mi/2) and /x > 0. We use mt = 178 GeV as suggested 
by the latest CDF and DO results [15]. We see that Bq generally increases monotonically 
for all positive values of Ao, and also for some negative values of Ao- This is also true for 
H < 0, as seen in Fig. 2. These observations immediately imply that, in any VCMSSM that 
predicts a unique value of Bo for a given value of Ao, there will be (at most) a unique value 
of tan /3 where the VCMSSM relation is obeyed. 
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Figure 1: Values of Bq as functions of tan j3 for n > and (mi/2,'n^o) = T^j (200, 1000) GeV, 
(b) (1000,1000) GeV, (c) (200,200) GeV and (d) (1000,200) GeV. Solutions for B = A- I 
case are denoted by small circles, which are connected by dashed lines. Solutions in the case 
of the Giudice-Masiero mechanism are denoted by small squares, connected by dot-dashed 
lines when possible. 
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Figure 2: As in Fig. 1, but for jj, < 0. 
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We note, however, that there are some particular negative values of Ao for which the 
required value of Bq, after rising when tan/3 is small, decreases slightly at large tan/3. This 
raises the possibility that there might be two allowed values of tan (3 in some restricted set 
of VCMSSMs. One example is when mi/2 = 200 GeV, mo = 1000 GeV, // > 0, i ~ -2 and 
B ~ —1.5, as seen in panel (a) of Fig. 1, and there are some other examples in other panels 
of Figs. 1 and 2. However, in practice, such multiple solutions do not exist in the specific 
VCMSSMs that we study in this paper. 

To illustrate this more explicitly, we indicate by small circles in Figs. 1 and 2 the values 
of tan (3 where the minimal supergravity condition B — A — 1 is satisfied, for a few specific 
values of A, and we join these points by dashed hues. For each value of A there is clearly 
only one consistent choice of tan /3 for any given choice of (mi/2, mo). We also show solutions 
for the Giudice-Masiero mechanism case, indicated by small squares. 

Note that we do not obtain solutions for B for all choices of A. For example, in Fig. la, 
we show solutions only for ^4 = 0, 1 and 2 ior B — A—1 and ^4 = 1 for the Giudice-Masiero 
model. In the minimal supergravity case, when A is reduced, B is also reduced driving 
the solution to smaller values of tan /3. Very quickly these solutions drop below tan P — 2 
and, below tan/3 ~ 1.7, the RGEs do not provide solutions to the sparticle spectra due to 
a divergence in the top quark Yukawa coupling at the unification scale. Similarly when A 
is large, the solution is driven to very large values of tan /3 where again no solutions to the 
RGEs are found. In the case of the GM model, the slope of Bq vs tan /3 is very small, and 
small changes in A lead to large changes in tan /3. Note also that in the GM model, there are 
often two branches of solutions which are disconnected. This is seen for example in Figs. Id 
and 2d. This is due to the relation (1) which separates solutions a,t A — 3. 

In order to have an analytical feel for the solutions for Bq shown in Fig. 1 and 2, we show 
in Fig. 3 the values of Bq at the input GUT scale, the tree-level values at the electroweak 
scale and the full values of B{Mw) as functions of tan /3, (a) for // > and (b) for // < 0, 
in both cases for (mi/2,mo) = (200,200) GeV. The tree-level value of B at the electroweak 
scale is defined as 

_ (mf + mi + 2/x^)sin2/3 

J^tree = '^^2^ ' 

and as one can see Btree tends to as tan/3 is increased. The 'full' values are calculated 

including one-loop electroweak threshold corrections, and Bq is then the result of running 

the RGEs from the weak scale to the unification scale. In the n > case, we see in Fig. 3 

that Bq is systematically larger than the tree-level value of B{Mw), which is in turn larger 

^As observed in [6], comparisons [16] with ISASUGRA [17] show that our procedure of minimizing the Higgs 
potential at the weak scale gives very similar spectra, also at large tan f3 and in the focus-point region. 
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Figure 3: The values of Bq, the electroweak tree-level and full B{Mw), as functions of tan f3 
for (a) ii> and (b) fi<0, both for (mi/2, mo) = (200, 200) GeV and A = l. 

than its full value. However, even in this case B{Mw) increases monotonically with tan/?. 
The situation is rather different for /i < 0, where we see that the sign of the loop correction 
depends on the value of tan (3, vanishing for tan /3 ~ 21. As a result, the full value of B{Mw) 
and hence Bq increase monotonically with tan f3. Had we neglected the 1-loop corrections to 
B and ran the RGEs up to the unification scale, we could obtain a non-monotonic solution 
for Bq with respect to tan (3 (for example, a solution with a minimum value of Bq) leading to 
multiple solutions of tan/5 for a fixed value of A [18]. Thus, Fig. 3 indicates the importance 
of the loop correction in determining the number of solutions. 

4 Updated Constraints on the CMSSM 

The standard LEP constraints and cosmological constraints on the CMSSM have been dis- 
cussed previously in many places [19, 20], so we do not discuss them further here, except to 
recall that we use the WMAP range 0.094 < ^l^h'^ < 0.129 [21] for the relic density of the 
LSP, assumed to be the lightest neutralino x- However, there are three new experimental 
developments that we should like to mention. One is the new value rrit = 178.0 ± 4.3 GeV 
recently reported by the CDF and DO collaborations [15], another is the evolution in the 
possible discrepancy between the experimental measurement oi — 2 and the value calcu- 
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lated in the SM, and the other is a recent improved upper hmit on the branching ratio for 

The new value of rrit affects the CMSSM parameter space in three important ways. One 
is to alter the calculation of the lightest MSSM Higgs boson mass rrih, and hence the lower 
limit on mi/2 inferred from the LEP lower limit ruh > 114.4 GeV. For example, in panel (a) 
of Fig. 4 for tan P — 10 and > 0, the lower limit on mi/2 is reduced by about 50 GeV when 
one increases rrit from 175 GeV to the value of 178 GeV shown here. A second effect is to 
alter the calculation of the rapid-annihilation funnels shown in panels (c) and (d) of Fig. 4 
for tan 13 — 35 and // < and for tan P — 50 and > 0, respectively. The sensitivity of these 
regions to rrit and large tan /3 was discussed earlier [22] in the context of the observabihty 
of the Higgs boson at hadron coUiders. Finally, the larger value of rrtt increases significantly 
the value of mo where the focus-point region may be found [23] . For example, for tan (3 — 10 
and > 0, we now find a focus-point region only for mo ^ 7 TeV for mi/2 ^ 250 GeV. We 
do not discuss focus points further in this paper. 

The BNL — 2 experiment recently announced a new determination using /i" and a 
final combined value using all their data [24]. Comparing with the SM calculations of 
Davier et al. [25], they quote a discrepancy of = — 2)/2 with the SM amounting to 

5a^ = (27 ± 10) X 10-1° (e+e" data) 

= (12 ± 9) X 10-1° (r data). (16) 

Another calculation of the SM contribution to {g^ — 2) using just the e+e" data [26] yielded 
a slightly larger discrepancy: 

Sa^ = (32 ± 10) X 10-1° (e+e- data) (17) 

There has subsequently been a new SM calculation of the hadronic vacuum polarization 
contribution by de Troconiz and Yndurain [27], who quote 

Sa^ = (27 ± 8) X 10-1° (e+e" data) 

= (19 ± 8) X 10-1° (r and e+e" data). (18) 

However, neither of these evaluations include the recent re-evaluation of the light-by-light 
contribution to by Melnikov and Vainshtein [28] , which decreases the discrepancy with the 
SM by about 4 x 10-i° compared with (18). Therefore, for the purposes of the subsequent 
discussion, we show contours corresponding to 

Sa^ = (15 ± 8) X 10-1°. (19) 
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We exhibit this constraint at the 2-a level, in which case its effect is essentially to exclude the 
option < but allow most of the (mi/2, fno) plane for > 0, apart from a region of small 
mi/2 and mo- However, we are well aware that the range (19) is open to question, particularly 
in view of the discrepancy between the estimates of the SM contribution based on e"^e~ and 
T data, and, to a lesser extent, the uncertainty in the light-by-light contribution. Therefore, 
we use (19) only as an indication, and by no means a rigid constraint on the parameter space 
of the CMSSM or any VCMSSM. In particular, we do not discard the option // < 0. 

Finally, we note that the CDF Collaboration have recently pubhshed an improved exper- 
imental upper hmit on the branching ratio for Bg — > H'^jJT [29], namely 5.8 x 10~^. Since the 
branching ratio for this decay oc tan^ /3 in the CMSSM, this constraint is potentially impor- 
tant at large tan/3. We find that this constraint is currently still 'covered' by the constraints 
from h — > S7, and g^j^ — 2, but this situation may change in the near future. 

In preparing the (mi/2, mo) planes in Fig. 4 and the subsequent figures, we have updated 
our code by making improvements that have impacts principally in the rapid-annihilation 
funnels and focus-point regions ^. Their effects are smaller than the other effects mentioned 
above. 

We show in Fig. 4 the (mi/2,mo) planes for a popular set of CMSSM cases, namely 
(a) tan/3 = 10, // > 0, (b) tan/3 = 10, // < 0, (c) tan/3 = 35,// < 0, and (d) tan/3 = 
50, // > 0, all for A = ^. In each panel, as well as the 'standard' experimental and 
cosmological constraints, we have indicated some representative contours of B by (blue) 
dashed lines. We see that B ^ in almost all the (mi/2,mo) planes exhibited, which is 
clearly incompatible with the minimal supergravity condition B — A — 1. This exemplifies 
the point that parameter choices allowed in the 'standard' CMSSM are often not allowed in 
favoured VCMSSMs. Specifically, the CMSSM cases shown in Fig. 4 could not be realized 
in minimal supergravity: one needs to choose smaller values of tan /3. A similar conclusion 
applies to the Giudice-Masiero model, which would require B — 1 ior the case A — 
considered here, although CM solutions are possible for // < if one discards the gn — 2 
constraint and chooses tan /3 somewhat above 10. 

Fig. 5 shows the corresponding (mi/2,mo) planes for the same choices of tan/3 and the 
sign of fjL, but for A — -1-0.75. In this case, the laboratory and cosmological constraints are 
not greatly different, even in the rapid-annihilation funnel regions ^. Now, however, there are 

^Specifically, we now include the full one-loop corrections to nib and rut instead of approximate expres- 
sions [30], and we correct a minor coding error. 

^We note in panels (c) and (d) the appearance of allowed bands above the x~ coannihilation strips, 
which are due to rapid fifi H annihilation. 

"^We note again the rapid fifi H annihilation strips in panels (c) and (d). 
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Figure 4: The {mi/2,mo) planes in the CMSSM for (a) tan/? = 10, /i > 0, (b) tan/? = 
10, yU < 0, (c) tan/3 = 35, yU < 0, and (d) tan/5 = 50, yU > 0, all for A = 0. In each 
panel, the region allowed by the WMAP cosmological constraint 0.094 < Q^h"^ < 0.129 has 
pale (turquoise) shading. The disallowed region where < has dark (red) shading. 
The regions excluded by b s'-f have medium (green) shading, and those in panels (a,d) 
that are disfavoured by — 2 at the 2-a level have very light (yellow) shading with a thin 
(black) border. The contours m^± = 104 GeV (m^ = 114 GeV) are shown as near-vertical 
black dashed (red dot-dashed) lines. In addition, we show several contours of B as (blue) 
dashed lines. There is no allowed point compatible with the minimal supergravity condition 
B = A — 1 or the Giudice-Masiero model in ^ese plots. 



some points where the minimal supergravity condition B — A — 1 is satisfied, as shown by 
the intersection of the B — —0.25 fine with the WMAP coannihilation region in Fig. 5a. For 
this choice of A, the Giudice-Masiero model is also satisfied at a hmited number of points, 
exemplified by the intersection of the B — 2/3 line with the WMAP funnel region in Fig. 5d. 

5 Examples of (mi/2,mo) Planes in VCMSSMs 

We now discuss the impacts of the above constraints on some specific VCMSSMs within 
the general framework of minimal supergravity, in which B — A — 1. As usual, we display 
these constraints in (mi/2, ttiq) planes. For the reasons discussed earlier, we regard tan/3 as 
a dependent quantity that varies across these planes, rather than being a fixed quantity as 
in most CMSSM analyses. Another difference from most CMSSM analyses is that the latter 
generally consider only the possibihty that the LSP is the hghtest neutralino x, assuming 
implicitly that the gravitino mass 1113/2 is sufficiently large that the gravitino LSP possibility 
can be neglected. However, in minimal supergravity, one has 1713/2 — mo (4) if the cosmo- 
logical constant A = 0, and the identity of the LSP varies over the (mi/2,mo) plane. We 
have recently published an analysis which includes the possibility that the gravitino is the 
LSP possibihty [8] , taking into account the constraints imposed by Big-Bang nucleosynthesis 
(BBN) and the cosmic microwave background (CMB) data on decays of the next-to-hghtest 
sparticle (NSP) into the gravitino, as well as the rehc gravitino dark matter density itself. 
In this paper, we incorporate this analysis into a unified treatment of the neutralino and 
gravitino LSP regions of the (mi/2,mo) planes in VCMSSMs. 

We display in Fig. 6 the contours of tan/3 (solid blue lines) in the (mi/2,mo) planes for 
selected values of A, B and the sign of Also shown are the contours where m^± > 104 GeV 
(near- vertical black dashed fines) and > 114 GeV (diagonal red dash-dotted lines). The 
regions excluded by 6 — > 57 have medium (green) shading, and those where the relic density 
of neutrafinos fies within the WMAP range 0.094 < Q,^h'^ < 0.129 have light (turquoise) 
shading. The gravitino LSP and the neutrafino LSP regions are separated by dark (chocolate) 
solid fines, and the WMAP relic-density strip for neutrafinos is shown only above these lines. 
The regions disfavoured by gi^u — 2 at the 2-cr level are very fight (yellow) shaded. 

If A has a large negative value, we do not find any consistent solutions to the electroweak 
vacuum conditions. This is refiected in panel (a) of Fig. 6, for > and A — —1.5, where 
there are no solutions above the topmost sohd (black) fine. The sofid (blue) contours of tan /3 
rise diagonally from low values of {mi/2,mo) to higher values, with higher values of tan/3 
having lower values of mo for a given value of mi/2. The dash-dotted (red) = 114 GeV 
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Figure 5: The (mi/2,mo) planes in the CMSSM for the same values of tan (3 and the sign of 
fi as in Fig. 4, but here for A = +0.75. 
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Figure 6: Examples of VCMSSM (mi/2,mo) planes with contours of tan P superposed, for 
/i > and (a) A = —1.5, (h) A = 0.75, (c) the simplest Polonyi model with A = 3 — ^/S, 
and (d) A = 2.0, all with B = A — 1. In each panel, we show the regions excluded by the 
LEP lower limits on MSSM particles and those ruled out by b ^ sj decay [31] (medium 
green shading): the regions disfavoured by — 2 are very light (yellow) shaded, bordered by 
a thin (black) line. The dark (chocolate) solid lines in panels (a, b, c) separate the gravitino 
LSP regions (below). Panel (d) exhibits a dark (red) wedge where the LSP is the rapid fi. 
The regions favoured by WMAP in the neutralino LSP case have light (turquoise) shading. 
The dashed (pink) line corresponds to the maximum relic density for the gravitino LSP, and 
regions allowed by BBN/CMB constraint on MSP decay are light (yellow) shaded. 



contour rises in a similar way, and regions above and to the left of this contour have ruh < 114 
GeV and are excluded. In particular, a neutralino LSP is excluded in this case. We exhibit 
in this and the other panels a gravitino LSP region, which was not studied in our previous 
exploration of VCMSSMs [6]. The relic density is acceptably low only below the dashed 
(pink) line. This excludes a supplementary domain of the (mi/2, mo) plane, but the strongest 
constraint is provided by the BBN/CMB decay constraint (hght, yellow shading), which 
requires tan/3 ^ 4.5. In panels (b, c, d), the rrih contour rises more vertically, but only in 
panel (d) is there any allowed neutralino LSP region. Panel (d) features an excluded dark 
(red) shaded wedge where the LSP is the fi. 

When A is increased to 0.75, as seen in panel (b) of Fig. 6, both the tan /3 and ruh contours 
rise more rapidly with mi/2- Again, there is no allowed neutralino LSP region. Within the 
gravitino LSP region, the rrih and relic density constraints would both be compatible with 
tan/3 ^ 7.5, but the BBN/CMB decay constraint imposes the stronger constraint that 
tan /3 ^ 13. It is instructive to compare this figure with Fig. 5a, which both assume that 
A — 0.75. The most notable difference is that, here, fixing the gravitino mass to equal mo 
excludes the neutralino coannihilation region with B — —0.25 and allows a region of the 
(mi/2,mo) plane that would previously have been excluded because the LSP would have 
been the fi. 

An analogous pattern is seen in the simplest Polonyi model with ^ = 3 — -\/3 shown in 
panel (c) of Fig. 6, where we note that the tan (3 contours have noticeable curvature. Once 
again, the neutralino LSP region is excluded, now by a combination of the Higgs and chargino 
mass bounds. At low mo in the gravitino LSP region, the rrih and relic gravitino density 
constraints impose tan/3 ^ 10 and the BBN/CMB decay constraint imposes tan/3 ^ 14 ^. 

We consider finally the case A — 2.0 shown in panel (d) of Fig. 6. In this case, there is a 
neutralino LSP region in the x — f coannihilation strip. Without the g'^u — 2 constraint, the 
most severe constraint on this region is imposed by 6 — > sj, requiring tan/3 ^ 25, the rrih 
constraint being much weaker. Imposing the fl'^u — 2 constraint requires tan /3 ^ 27. There is 
an excluded dark (red) shaded wedge where the LSP is the fi. Below this appears a gravitino 
LSP region with acceptable rehc density. Within this region, the rrih and BBN/CMB decay 
constraints impose tan/3 ^ 15, which would be strengthened to tan/3 ^ 20 if one took the 

— 2 constraint at face value. This is the shaded region in the lower right of panel (d) . 

We find no consistent solutions for values of A substantially greater than 3 (4) when 
> (// < 0), and negative values of A are not allowed when fi < 0. These restrictions arise 

^There is also a negative Polonyi solution with A = —3 + \/3, whose {mi/2,mo) plane is qualitatively 
similar to panel (a) of Fig. 6. 
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Figure 7: As in Fig. 6, but now for /i < and the choices (a) A = 3 — B = A — 1 and 
(b) A = 2,B = A-1 and fi<0. 

from the behavior of the relation between tan/5 and Bq discussed earher. Therefore as A 
increases, so does the solution for tan/5 when B = A — 1. At very large tan/5, there are no 
solutions to the RGEs due to a divergence in the bottom-quark Yukawa coupling. For small 
A and /i < 0, the solution is driven to excessively small values of tan/5, where again there 
are no solutions, now due to the divergence in the top Yukawa coupling mentioned earlier. 
The same is true when A is large and negative and /i > 0, i.e. for A < —2.5, mo ^ 500 GeV 
for mi/2 < 1000 GeV. 

As we see in panel (a) of Fig. 7, only a small area of the {mi/2, mo) plane in the gravitino 
LSP region is allowed by the mh constraint in the positive Polonyi case i = 3 - 73 6 This 
area would be further restricted if one took the Qfi — 2 constraint at face value. At larger 
values of A, the allowed region is extended, as exemplified in panel (b) of Fig. 7 for the 
case A = 2, where the rrih constraint is somewhat weaker. However, in this case the Qfj, — 2 
constraint would have a much more drastic effect. 
^The negative Polonyi case is not allowed for /i < 0. 
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6 The Problem and the Giudice-Masiero Mechanism 



One of the primary motivations in building supersymmetric model is to avoid the the gauge 
hierarchy problem, namely that the Higgs mass is of order mz and much less than the Planck 
mass, though not protected by any symmetry of the Standard Model between the GUT scale 
and the weak scale. Supersymmetry alleviates this problem via cancellations between contri- 
butions to the Higgs mass from fermions and bosons in the same supermultiplet. However, 
this scenario begs the question why supersymmetry is broken by soft terms which are as- 
sumed to be 0(1 TeV). Moreover, there is one other, supersymmetric, parameter which is 
required to be small, namely the Higgs mixing parameter /x. One of the most interesting 
attempts to explain the smallness of ji is the Giudice-Masiero mechanism [5], in which it 
is related to a coupling between observable and hidden sectors, and is of the same order 
of magnitude as the soft supersymmetry-breaking parameters. In the simplest realization 
of the Giudice-Masiero mechanism with only one hidden superfield, one has the following 
relation between A and as already mentioned: 

B = H^, (20) 



and 



mo 



V3 



(21) 



where A is the coupling constant between the hidden superfield and the two Higgs supermul- 
tiplets. One should require that A ~ 0(1) for fj, to be the same order of mo- 

We display in Fig. 8 some typical (mi/2,mo) planes in the Giudice-Masiero model for 
positive fjL. As in the previous minimal supergravity VCMSSMs, we find no consistent elec- 
troweak solutions for values of A much outside the range of values exhibited. In the examples 
shown, there are no solutions above the topmost solid lines in panels (a) for A — 0.6 and (d) 
for A = 1.8. For A < -0.6, mo ^ 150 GeV for mi/2 < 1000 GeV. Similarly for A ^ 2.6 only 
a small corner of the plane admits solutions. 

In panel (a) for A — 0.6, corresponding to B — 0.75, there is no allowed region above the 
solid (chocolate) gravitino LSP hue. Below this line, we see an allowed region for tan /3 ^ 22. 
However, we also note that the corresponding values of A are quite large, A ^ 5. The 
situation is somewhat different for the case A — 0.8, corresponding to B ~ 0.64, shown in 
panel (b) of Fig. 8. In this case, we see that there is a narrow allowed region along the 
X — f coannihilation strip in the neutralino LSP region for tan/3 ^ 33.5, or ^ 35 if the 
Qix — 2 constraint is taken into account. This region requires A ^ 2, which is relatively 
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palatable. At lower mo, there is a disallowed dark (red) wedge where the fi is the LSP, and 
below that a region where the gravitino is the LSP. The latter contains a domain allowed 
by the BBN/CMB decay constraint, that appears for tan/3 ^ 18, or ^ 20 if one includes 
the — 2 constraint. However, this region again has A ^ 5. Turning now to the Polony i 
case A — ?, — -\/3 shown in panel (c) of Fig. 8, corresponding to S ~ 0.27, we see that there 
is no allowed area in the neutralino LSP region above the dark sohd (chocolate) line, and 
that the allowed region in the gravitino LSP region requires tan/3 ^ 13.5 and again A ^ 5. 
Similar features are seen in panel (d) for A — 1.8, corresponding to B — —0.5, where the 
only allowed area - in the gravitino LSP region - requires even larger values of A than the 
previous cases. 

Fig. 9 shows some analogous cases for < 0. As before, there are no consistent elec- 
troweak vacuum solutions for values of A substantially outside the range of values shown, and 
none above the topmost solid lines in panels (a) and (b). Panel (a) is for A — —0.2, corre- 
sponding to S ~ 1.06. It has two narrow strips in the neutralino LSP region that are allowed 
if one discards the 9^ — 2 constraint, appearing for mi/2 ^ 800 GeV and mo ^ 500 GeV 
for tan/3 > 38 and A < 1.5. Down in the gravitino LSP area, there is a second allowed 
region with tan /3 ^ 13 and A ^ 3. For smaller values of A, the allowed parameter space is 
further squeezed. For example, for A — —1 we find mo ^ 500 GeV for mi/2 < 1000 GeV. 
In panel (b) for A — 0.6, corresponding to B — 0.75, the allowed neutrahno LSP region has 
disappeared, but a gravitino LSP region remains. Similar features are seen in panels (c) and 
(d) for i = 1 (S = 0.5) and i = 3 - ^3 (S ~ 0.27), respectively. For i ^ 2, solutions exist 
only in a small portion of the plane. 

7 Conclusions 

We have discussed in this paper the impacts of the theoretical, experimental and cosmological 
constraints on some classes of VCMSSMs, including minimal supergravity models and the 
Giudice-Masiero model. We have presented unified treatments of the regions of parameter 
space in these models where the LSP is a neutralino or the gravitino. 

We have emphasized that the predictions of these models differ significantly from those 
of the CMSSM. In particular, the CMSSM is distinct from minimal supergravity: the for- 
mer does not necessarily require a fixed relation between the trilinear and bilinear soft 
supersymmetry- breaking parameters A,B, nor equality between mo and mz/2, as required 
in minimal supergravity models. The values of B required in generic realizations of the 
CMSSM generally bear no relation to the values that would be derived in minimal super- 
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Figure 8: Examples of {mi/2,mQ) planes in the Giudice-Masiero model, with contours of 
tan (3 superposed, for fj. > and (a) A = 0.6, (h) A = 0.8, (c) A = 3 — -\/3 and (d) A = 1.8. 
In each panel, we show the regions excluded by the LEP lower limits on MSSM particles, 
those ruled out by b —>■ S'j decay [31] (medium, green shading), and those disfavoured by 
— 2 (very light, yellow shading). As before, a dotted (red) line shows where = rrir and 
the gravitino LSP region is bounded by a solid (chocolate) line in panels (a, c, d). The dark 
(red) wedge in panel (b) has a fi LSP and hence is disallowed. Within the gravitino LSP 
region, the relic density constraint is indicated by a dashed (pink) line and the BBN/CMB 
constraint on NSP decay by light (yellow) shading. The dashed black lines are contours of 
the Giudice-Masiero parameter A. „„ 
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Figure 9: Same as Fig. 8, hut for /i < and (a) A = —0.2, (b) A = 0.6, (c) A= 1, and (d) 
A = 3-^/3. 



21 



gravity models. 

In addition to minimal supergravity models, we have discussed the simplest variant of the 
Giudice-Masiero model, which makes a brave attempt to provide a framework for calculating 
the Higgs-mixing superpotential term /i. 

There are a couple of striking features of these specific analyses that we note. One is that 
the range of A is often very restricted: beyond this range, it is impossible to find consistent 
solutions to the electroweak vacuum conditions. 

A second observation is that, in both minimal supergravity and the Giudice-Masiero 
model, a neutralino LSP is completely excluded in many instances, and the gravitino LSP 
regions are generally much more extensive than the neutralino LSP regions. To some extent, 
this was to be expected, since we impose the cosmological dark matter density and NSP decay 
constraints on gravitino dark matter as one-sided upper limits, rather than as narrow WMAP 
ranges as for the dark matter density constraint on neutralino dark matter. This is because, 
in the case of gravitino dark matter, the narrow range could be reached by postulating 
thermal gravitino production with a suitable reheating temperature [32]. Of course, in 
either the neutralino or gravitino case, one could always postulate a supplementary source of 
cold dark matter. In the case of neutralino dark matter, this possibility would broaden the 
WMAP density strip down to the — ruf boundary. However, the gravitino dark matter 
region would still, for many choices of the other supersymmetric parameters, occupy a larger 
area of the (mi/2, mo) plane. 

In any complete supersymmetric theory, one expects some relations between supersym- 
metry breaking parameters, perhaps of the type discussed here. In this case, some VCMSSM 
should be responsible for the low-energy sparticle spectrum. However, we do not yet know 
what specific constraints are handed down from the unification or string scales. As we have 
emphasized in this paper, the predictions in such models may differ greatly from those of the 
more relaxed CMSSM and, a priori, those of a more general MSSM. Analogous differences 
are also to be expected in the predicted cross sections for direct and indirect searches for 
supersymmetric dark matter, a topic we will consider elsewhere. 
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